Nature of the inhomogeneous state of the extended t-J model on a square lattice 
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We carry out the variational Monte Carlo calculation to examine spatially inhomogeneous states in 
hole- and electron-doped cuprates. By using Gutzwiller approximation, we consider the excitations, 
arising from charge density, spin density and pair field, of the mean-field ground state of the t — J 
model. It leads to the stripe patterns we have found numerically in a generalized t — J— type 
model including mass renormalization from the electron-phonon coupling. In the hole-doped side, 
a robust d-wave superconducting order results in the formation of the half-doped antiferromagnetic 
resonating- valence-bond (AF-RVB) stripes shown by the well-known Yamada plot. On the other 
hand, due to a long-range AF order in electron-doped materials, a stripe structure with the "in- 
phase" magnetic domain (IPMD) is obtained in the underdoped regime instead of the AF-RVB 
stripe. The IPMD stripe with the largest period permitted by lattice size is stabilized near the 
underdoped region and it excludes the Yamada plot from electron-doped cases. Based on finite 
lattice size to which we can reach, the existence of IPMD stripes may imply an electronic phase 
separation into an electron-rich and an insulating half-filled AF long-range ordered domains in 
electron-doped compounds. 

PACS numbers: 74.20.-z, 74.72. Ek,71.10.-w,71. 10. Fd 



Since the discovery of high-temperature superconduc- 
tivity in the layered cuprate materials, there have been 
many evidences for stripe structures in several families of 
the hole-doped cuprates, e.g., La2-sSrsCuOr^. One of 
the many puzzles of the stripes is the doping dependence 
of the incommensurate magnetic peaks associated with 
the stripes measured by neutron-scattering experiments 
which obeys the so-called Yamada plot-. It represents 
the existence of the half-doped stripe with average of 1/2 
hole in one charge modulation period at 1/8 hole density 
and below. There have been many early theoretical works 
attempting to explain the Yamada plof*^^. One possible 
scenario for such correlation is that tendency for charges 
toward phase separation can lead to various structures 
including stripes, puddles^i^, or even cluster glasses with 
randomly-oriented stripe domains recently observed by 
scanning tunneling spectroscopy (STS)^'^°. Recently we 
have used a variational Monte Carlo (VMC) technique^i 
to successfully establish the half-doped stripes in the ex- 
tended t — J— type Hamiltonian by including a mass- 
renormalization effect due to a weak electron-phonon 
coupling. 

So far the experimental situation in electron-doped ma- 
terials is much less cleariS,. There are several indirect ev- 
idences for a homogeneous state in electron-doped com- 
pounds coming from measurements of neutron scatter- 
ing and core-level photoemission spectraii^iii^. Yamada 
et al. have reported only commensurate spin fluctua- 
tions observed by neutron scattering in the superconduct- 
ing (SC) Ndi,s5Ceo.i5CuO^. It is different from the 
incommensurate peaks observed in hole-doped cuprates 
which are considered as the hall mark of the "out-of- 
phase" stripe domains with a 7r-phase-shifted staggered 
magnetic moment. Instead of stripes, the short-range 
spatial inhomogeneity of the antiferromagnetic (AF) cor- 
relations was recently reported for the electron-doped su- 
perconductor Pro,gsLaCeo,i2Cu04-s by Zhao et al^*^ by 



using STS and neutron scattering. In addition, there 
are also evidences to support inhomogeneous states from 
measurements of muon spin rotation (/iSR}^'^^, nuclear 
magnetic resonancei^, magnetoresistance^S, and ther- 
mal conductivity^^. Whether Inhomogeneity in electron- 
doped compounds is intrinsic or induced by the cerium 
doping or oxygen defects were discussed by two recent 
works^^'^'^. All these results suggest that the possibility 
of phase separation and inhomogeneity is an unresolved 
issue in electron-doped cuprates. 

On the theoretical side, there are some numerical ev- 
idences that Hubbard models produce stripes in the 
electron-doped system. Within an unrestricted Hartree- 
Fock approach, the occurrence of the diagonal filled 
stripes having average of one doped electron per stripe 
site has been demonstrated earlier—. Later, the vertical 
"in-phase" stripe domains without 7r-phase-shifted stag- 
gered magnetic moment in the t — t' Hubbard model 
have been found in the electron-doped regime^. The 
unusual doping evolution of the Fermi surface detected 
by angle-resolved photoemission spectroscopy^^ was ex- 
plained by assuming the inhomogeneous in-phase stripe 
phases?'' but without including strong correlations. Since 
the strong correlation makes the difference in energies 
between uniform states and various hole-doped stripe 
states very small which has been shown recently^^''^^, it 
is beyond the accuracy of Hartree-Fock method to ad- 
dress this kind of difference. Therefore, a careful varia- 
tional approach is needed to examine the stability of the 
electron-doped stripe states. 

In this paper, we shall first demonstrate that the par- 
ticular spatial patterns of modulations of the charge den- 
sity, spin density and pair field could be derived by con- 
sidering the excitations of the mean-field ground state of 
the extended t — J model using Gutzwiller approxima- 
tion. Once the relations between charge, spin and pair 
field are revealed we then explicitly construct the stripe 
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wave functions. As a comparison with what we have done 
previously for the hole-doped cases^, we shall add an 
electron-phonon interaction to the model before carrying 
out the numerical calculations. Just as before, we will 
not consider the full effect of electron-phonon coupling 
but only examine the simplest effect of mass renormal- 
ization of charges due to phonon couplings. The renor- 
malization effect depending on the local carrier density 
is treated self-consistently in the VMC method taking 
into account the strong correlation exactljJ^. We find 
that half-doped stripes obtained in the hole-doped sys- 
tems are no longer stable in the electron-doped cases for 
a range of electron-phonon interaction strength. Instead 
the system is very likely to have an electronic phase sep- 
aration. The different results between hole-doped and 
electron-doped systems is mainly due to the strong AF 
long-range order in the latter system. The effect of t' /t 
will be also discussed. 

We consider the extended t — J Hamiltonian on a 
square lattice given by 



H = 



+ H.C. 



J ^ Si • , 
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where the hopping tij — t, t' , and t" for sites i and j being 
the nearest, second-nearest, and third-nearest neighbors, 
respectively. Other notations are standard. In the follow- 
ing, the bare parameters t" and J in the Hamiltonian are 
set to be {t",J)/t = (— 1'/2, 0.3). Since doubly-occupied 
sites in electron-doped materials play the same role as 
holes in hole-doped cases, we treat the hole- and electron- 
doped cases in the same manner except that t' /t — ?> —t'/t 
and t"/t -t"/i^. t'/t < (> 0) corresponds to the 
hole-doped (electron-doped) regions. In this paper, we 
primarily study the hole- and electron-doped phase dia- 
grams for different t'/t. 

In a generalized mean-field theory to include the pos- 
sibility of spatially non-uniform solutions, we define the 
local carrier density pi, the local AF order to^, and the 
nearest- neighbor pair field . The mean- field Hamilto- 
nian is simply given by 



MF 
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Si 
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where the matrix elements 



A, = 



(3) 



Here 1, 2, and 3 correspond to the nearest, second- 
nearest, and third-nearest neighbors, respectively, and 
a=t(l) ori(-l). 

Once the variational parameters p^, nii, and Aij are 
given, we can diagonalize Eq. ([2]) to obtain N positive and 
iV negative eigenvalues with corresponding eigenvectors 



(Mf,wf) and (u",t;"), given by 



(4) 



Here A'' is the lattice size. We can formulate the trial 
wave function fixing the number of electrons with 
the Gutzwiller projector Pq and the hole-hole repulsive 
Jastrow factor Pj (sec the details of Ref . [28l) . 

IvP) = PgPjPn^HiuiHo) 
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To avoid the divergent determinant because of the 
presence of nodes in the RVB-type wave functions 
with periodic boundary condition, a particle-hole 
transformation^^iM, c-^ fi and c-^ -> dj, has been 

introduced in Eq.®. Here \Q) =IU,\0). 

In principle, pi, m^, and could be determined vari- 
ationally. In practice, it is very difficult to optimize the 
energy of such multi-variable problems. If we postulate 
that the inhomogeneous states are actually fluctuations 
beyond the uniform mean-field solution, then a more effi- 
cient method to find the most probable solutions is to use 
Gutzwiller approximatioi*^. In this approximation, the 
total energy (H) with respect to the Gutzwiller-projected 
wave function can be written as 



E U,gn^)9nJ) {X^j + H.c.) 
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where the Gutzwiller factors gt{j) and gs{i) are known 
to be2i 



9s{i) 
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Here = 

X*j(= Y.aXfj = E<r(4<TCjT)o), {Sf)o), and Ay(= 

{ciiCjf — Ci^Cj i)o) is the bond order parameter, staggered 
magnetization, and pair field with respect to the non- 
projected wave function, respectively. For a usual mean- 
field theory, these parameters are assumed to be constant 
and the same for all sites or bonds. Here we shall go one 
step further by examining the fluctuations beyond the 
constant mean-field values as given by 



6. 
m. 
A 



S + dSi 
ffi + dm. 
A + dA. 



(8) 
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Here d{...) means the small fluctuation away from the 
average value (...). By substituting Eq.® into Eq.®, 
we obtain several coupled terms with the form dSidmf, 

There are also non- 



and dSid Af^ 
coupled terms of the form dm| 



AdSidAij, mdSidm. 



ddiddj, and dAf^. We 



neglect the fluctuation of bond order and also skip the 
derivation of all the terms as they are irrelevant for our 
calculations below. 

In the hole-doped side, at flnite doping there is no long- 
range AF order (m = 0) but with a nonzero d-wave SC 
order parameter (A ^ 0). Thus, according to the discus- 
sion above, the relevant fluctuation terms to couple the 
charge, spin and pair flelds are dSidmf and AdSidAij. 
The Fourier transform of these two terms are dSadm? 



and AdSqdA. 
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q, respectively. Hence the modulation pe- 



riod of charge density, Oc, should equal to the period of 
pair field, ap, but is only half the period of the spin den- 
sity, a^, i.e., ac = ap = as/2. More precisely the wave 
functions to include these fluctuations can have order pa- 
rameters of the form: 



Pv cos 



27r 



M ■ 

m„ sm 



A^,i+£ = A„ cos 
A, 









\- 

as 






"27r 






ap 




^)1 



A 



A5 cos 



2-K 



-Vi 



(9) 



Here ac — ap — as 12. The variational parameters are 
indicated by subscript "w". The average charge den- 
sity is determined by the chemical potential not in- 
cluded in Pi. This is exactly the wave function called 
the AF resonating-valence-bond (AF-RVB) stripe state 
by us in Ref. HI- There we had shown that the varia- 
tional energy of the uniform d-wave RVB (d-RVB) state 
(p„ = TO*^ = A*-'^ — 0) is considered as the reference 
energy. Furthermore if we add a weak electron-phonon 
interaction to the extended t — J model, as shown in 
Ref. [Ill, the half-doped AF-RVB stripes have lower en- 
ergy than the d-RVB state. This AF-RVB stripe pattern 
is in good agreement with experiments as well^'^^. 

The electron-phonon interaction is introduced by as- 
suming the hopping terms tij in Eq.([T]) modified due to 
the spatial variation in carrier density in the sense that 
the sites with larger carrier density in the modulated 
charge pattern have larger Uf^. Then we assume a linear 
relation between the electron-phonon coupling strength 
A and doping density S, X{S)/X{0) = 1-36 = f{S). Now 
tj^ can be renormalized to 
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where is the carrier density at site i and A the bare 
parameter in the Hamiltonian. The details are given in 
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FIG. 1: (Color online) The doping dependence of the op- 
timized energy difference in terms of percentage for several 
trial states denoted in the figure. The reference ground state 
is the uniform d-RVB state. The positive (negative) doping 
density S means the hole-doped (electron-doped) cases. Note 
that the AF-RVB stripe states with different magnetic peri- 
ods are indicated. The lattice size V is 24 x 24. Here A = 0.25 
and \t'/t\ = 0.1. 



Ref. [Til In what follows, we will use a notation "t* — J" 
to stand for this Hamiltonian. 

In the electron-doped side, it had been shown^ that 
there is a robust long-range AF order (m ^ 0) in the 
underdoped regime in contrast to hole-doped systems. 
Based on our previous discussion of fluctuations using the 
Gutzwiller approximation, the important contributions 
in Eq.® will be rhdSidrrii and AdSidAij. Hence, the 
magnetic modulation has the same period as charge and 
pair field, that is, ac = ap = as. This pattern is very 
different from AF-RVB stripe states discussed above for 
hole-doped systems. The wave function to describe such 
a new pattern is called the in-phase magnetic domains 
(IPMD). The IPMD stripe state has the same function 
of Pi and Aij as Eq. © , only the magnetic modulation is 
now changed to 



-mt^ cos 



2n 



(11) 



where a finite staggered magnetic moment m'^ is in- 
cluded. Due to the presence of long-range AF order for 
electron-doped systems, we have four variational states 
to be considered. These are the pure d-RVB uniform SC 
state, the AF-RVB stripe state, the IPMD stripe state 
and the uniform state with the coexistence of long-range 
AF and SC orders (AF-hSC)^. The AF-hSC state can 



in 



be simply constructed by setting p^ 
Eqs.® and (fTT|) . 

Figure [T] shows the percentage of energy change with 
respect to the d-RVB state for three other trial wave 
functions in the extended t* — J model with smaller 
\t' /t\{= 0.1) for hole- and electron-doped phase diagrams. 
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FIG. 2: (Color online) The profile of carrier density ns and 
staggered magnetization (M) for the optimized states: (a) the 
AF-RVB stripe and (b) the IPMD stripe with 24ao magnetic 
periodicity at 1/12 electron-doping in the extended t* — J 
Hamiltonian. The bare parameters A = 0.25 and \t' /t\ = 0.1. 
All quantities are calculated in a 24 x 24 lattice. 

In the hole-doped case, the half-doped AF-RVB stripe is 
still observed like our previous studies for t' /t{— — 0.2)-i. 
In fact, the half-doped AF-RVB stripe can be always 
found for all \t' /t\ we have investigated in the hole-doped 
region as shown in Fig|31[a) and (b). In other words, the 
Fermi surface topology seems not to affect the stability 
of the half-doped stripe obtained in the extended t* — J 
model. 

In the left panel of FigUwe find the IPMD stripe states 
are stabilized for electron- doping less than 0.1. For dop- 
ing greater than 0.1, the AF-RVB stripe state with the 
largest magnetic period has the lowest energy (24ao is 
the largest size we have studied here). The magnetic pe- 
riod does not change with the doping density as the half- 
doped stripes. We also examine the IPMD stripe states 
with different periods in addition to 24ao. The results 
are not shown here, but the most stable IPMD stripe 
has the largest magnetic period as lattice size. The dif- 
ference between the AF-RVB and IPMD stripes for the 
maximum magnetic period of 24ao is actually negligible. 
Since the IPMD state has much lower energy for doping 
less than 0.1, it appears that the t* — J Hamiltonian does 
not prefer to break the system into many 7r-phase-shift 
magnetic domains. In Figl^a) and (b) for 1/12 elec- 
tron density, the spatial variation of charge density and 



FIG. 3: (Color online) The same descriptions are presented 
in the caption of Fig[l] except for (a) \t' /t\ — 0.2 and (b) 
\t'/t\ = 0.3. 



staggered magnetization along the direction of modula- 
tion are shown for the AF-RVB stripe and IPMD stripe, 
respectively. For both states, there is almost no doped 
carriers in the region with the strongest staggered mag- 
netization. It is essentially a phase separated state with 
an electron-rich region and an insulating half-filled AF 
long-range ordered region. The IPMD stripe state has 
staggered magnetization |(Af)| ~ 0.265, which is close to 
the homogeneous AF-hSC state (|(M)| = 0.269). This is 
consistent with the commensurate short-range spin cor- 
relations detected by neutron scatteringi^. 

For larger t'/t, we expect robust AF correlations 
should persist to the higher doping in electron-doped 
cases. Figure [3] shows this common feature for the uni- 
form AF+SC and the IPMD stripe states. In Figl^^a), 
while the homogeneous AF-I-SC state has much lower en- 
ergy than the d-RVB state, the IPMD stripe state is 
still the best candidate for electron-doping (5 < 0.1. It 
is worth pointing out that for electron-doping 5 > 0.1 
the IPMD stripe state begins to gain energy to compete 
with the AF-RVB stripe state in the case of t'/t ~ 0.2. 
Interestingly, as increasing t'/t{— 0.3) further, the AF- 
RVB stripe states entirely disappear in the electron- 
doped phase diagram possibly due to its 7r-phase-shift 
domains in the spin modulation, as shown in FigjHJb). 
Except for the higher electron-doping region (S > 0.15) 
where the homogeneous AF-I-SC state still exists, the 
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FIG. 4: (Color online) A-dependence of the optimized energy 
difference in terms of percentage for several trial states de- 
noted in the figure at 1/12 doping. Here A'^ = 24 x 24 and 
t'/t = 0.3. 

IPMD stripe state is dominant within the wider doping 
range (0 < i5 < 0.15). In addition to these stripe states 
shown in the phase diagram, we have also examined other 
trial wave functions with different shapes and periods 
for stripes, such as glassy stripes with randomly-oriented 
8x8 magnetic patches, IPMD stripes with other periods, 
and the diagonal stripe with the largest magnetic period 
(12-\/2ao)- However, none of them can be stabilized for 
all t'/t in electron-doped cases (not shown). 

To complete the discussion on the VMC results, we 
study the effect of the electron-phonon coupling strength 
A on the phase diagrams. Since the IPMD stripe state 
has much lower energy near 1/12 electron-doping, we 
only show A-dependence of the optimized energy dif- 
ference for several trial wave functions including IPMD 
stripes with different magnetic periods at 1/12 doping 
in FigHl As we expected, the stability of the uniform 
AF-I-SC state is almost independent of A as A ^ 0.2. If 
the electron-phonon coupling strength is smaller (A < 
0.125), the AF-I-SC state would be still a good candidate 
for the ground state near the underdoped region in the 
electron-doped side. However, once A becomes larger, 
IPMD stripe states begin to be stabilized. Notably, the 
best one for A = 0.25 is either the IPMD stripe state 
with the largest magnetic period or the phase-separated 
state. Unlike hole-doped cases, there is no place for AF- 
RVB stripe states at this doping in electron-doped phase 
diagrams. 

Although the lowest energy solutions of electron-doped 
systems for the extended t* — J Hamiltonian are IPMD 
stripes and not AF-RVB stripes as the hole-doped sys- 
tems, the particular patterns of modulation periods of 
charge density, spin density and pair field for both stripes 
are derived from the fluctuations of the order parame- 
ters used in the mean field theory as shown by using the 



Gutzwiller approximation. The inhomogeneous states 
like stripes or phase separation is really due to excitations 
or fluctuations of order parameters of the ground states 
of the strongly correlated t — J model. When electron- 
phonon interaction is introduced to renormalize the mass 
of carriers, these inhomogeneous solutions then become 
stable. 

In conclusion, following the same approach we have 
used to study stripes for hole-doped systemaiii^, we have 
further studied the possibilities of having non-uniform 
ground states for the electron-doped cases using the VMC 
method. By using the Gutzwiller approximation to ex- 
amine the fluctuations of order parameters used to obtain 
the d-RVB ground states, we have found that the period 
of modulation of charge density, staggered moment and 
pair field are correlated. For hole-doped systems the low- 
est energy stripes, the AF-RVB stripe states, have the 
period of staggered moment twice of that of charge and 
pair field. These stripes with half a carrier per charge do- 
main become the ground states after we include the weak 
electron-phonon coupling. However, similar approach for 
electron-doped cases has turned out a different stripe pat- 
tern. The lowest energy stripes, the IPMD stripe states, 
now have the period of staggered moment same as that of 
charge and pair field. One of the main reasons is that for 
a large electron doping range the ground state has long- 
range AF order. The numerical results show that at low 
doping, the IPMD stripe state with the largest magnetic 
period same as lattice size is composed of an electron-rich 
region and an insulating half-filled long-range AF ordered 
region. We believe that it may indicate a phase separated 
state in electron-doped compounds. This statement re- 
quires the numerical calculation with large lattice size 
that is left for our future study. 

Our result that the electron-doped systems should not 
have the half-doped stripes observed in hole-doped sys- 
tems is consistent with the neutron scattering results for 
cuprates that Yamada plot is only observed for hole- 
doped cuprates^ but not for electron-doped^^. The evi- 
dence on magnetic inhomogeneity we have found numer- 
ically has been indirectly observed in the electron-doped 
cuprates Pr2-sCesCuOii and Pro,s8LaCeo,i2CuOi-s by 
using /iSR and STS experiments, respectively^^' The 
fact that for electron doping we predict possible phase 
separation which is also a hotly debated issue by experi- 
ments on cuprates is quite interesting. We look forward 
to possible resolution of this issue in the near future. 
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